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A full wave theory is applied to the whistler radio wave propagation along a plasma slab with an enhanced or depressed 
plasma density which is imbedded in an infinite magnetoplasma. Rigorous dispersion equation is solved for a thin slab 
in approximate but explicit forms. Three types of propagation modes are found: (a) a completely trapped surface wave 
mode along the depressed slab in the frequency region above a half of gyrofrequency and below a gyrofrequency, (b) a com- 
pletely trapped surface wave mode along the highly enhanced slab in the frequency region above a half of gyrofrequency 
and below a certain cutoff frequency less than a gyrofrequency, and (c) a partially trapped (leaky) surface wave mode 
along the enhanced slab in the frequency region above a certain cutoff frequency and below a half of gyrofrequency. Dis- 
persion properties, field distributions and an attenuation of the third mode due to the leakage of the transmitted power are 
discussed in detail. The attenuation is found to increase very rapidly with increasing frequency, thickness and enhancement 
of ionization of the slab. The exact numerical solutions are also obtained and compared with the approximate solutions. 



1. Introduction 

A number of theoretical investigations have been made on the propagation of whistler radio 
waves since the publication of Storey's famous paper [Storey, 1953], and not a few papers have 
dealt with the propagation path and the guiding mechanism of the whistlers. The present general 
idea accepted by many investigators on the guiding mechanism of the whistlers is that there may 
exist earth magnetic field-ahned ducts in which the electron density is greater (or less) than the 
density outside, and the whistlers are guided along the ducts without loss due to the divergence 
of the wave. Concerning the guidance along such ducts, good review papers have been given 
by Storey [1962] and Booker [1963]. 

The ray theory has been used for the calculations of the whistler ray paths by Maeda and 
Kimura [1956], Smith et al. [1960] and Yabroff [1961]. Voge [1961, 1962] and Booker [1962] 
calculated the minimum discontinuity of ionization density and the corresponding minimum 
width of the ducts required for the trapping by using the analogous method with that used in the 
tropospheric duct propagation. These theories are considered to be appropriate for compara- 
tively wide ducts. Northover [1959] has dealt with the guidance of the whistlers along an in- 
homogeneous cylindrical plasma column by use of the wave theory, in which he has expanded the 
fields into power series of the radial distance from the column axis. 

Storey [1962] has raised a question whether very narrow columns or slabs can support trapped 
waves analogous to the surface wave modes along open transmission hues. In order to answer 
this question a fuU wave theory must be applied. In connection with this question, this paper 
deals with the electromagnetic wave propagation along a sharply bounded homogeneous infinite 
magnetoplasma slab which is imbedded in an infinite magnetoplasma of a different plasma density 
with that of the slab. A rigorous dispersion equation based on the classical magneto-ionic theory, 
Maxwell's equations and conventional boundary conditions is derived and solved first for a thin 
slab in approximate but explicit forms. These approximate solutions are partly compared with 
numerically obtained exact solutions. The existence of completely trapped and partially trapped 
waves, their dispersion properties, and an attenuation due to the leakage are discussed in detail, 
in connection with the guiding mechanism of the whistlers. 
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Figure 1. A magnetoplasma slab imbedded in an 
infinite magnetoplasma of a different plasma density 
with that of a slab. 



2. Model for Theoretical Analysis 

The earth magnetic field-ahned plasma duct is represented by an infinite uniform plasma slab 
which is perpendicular to the jc-axis and has a thickness of 2a as shown in figure 1. The magnetic 
field of the earth is assumed to be uniform in the whole region and parallel to the slab and z-axis. 
The corresponding electron gyrofrequency of the magnetic field is denoted by o)//. Let the elec- 
tron plasma frequencies of the plasma slab and the surrounding plasma region be denoted by 0}p 
and ct)^, respectively. In the present analysis the effect of ions is neglected for simplicity. The 
effect of finite compressibility of the plasma is also disregarded. This does not cause any error 
when whistler radio waves propagating parallel to the magnetic field are considered. The error 
when the propagating direction is not parallel to the field is not appreciable unless the electron 
temperature exceeds lO'* °K [Kieburtz, 1964]. Thus the classical magneto-ionic theory is applied 
in this analysis. 

3. Electromagnetic Fields and Dispersion Equation 

The combined Maxwell's equation for an electric field E{x^ y, z) is expressed by 

(VV-V4-A:2|)£^=0, (1) 

where the time^ dependence of e^^^ is used and will be suppressed in the text, k^ is a free space 
wave number, 1 is the unit matrix, and e is the well-known relative dielectric tensor of the plasma 
and is given by the form 



J /e, -i€2 0^ 
e = jf.2 ei 

\0 e,,y 



(2) 



When the plasmas are assumed to be lossless, ei, €2, and 63 (the plasma parameters inside the 
plasma slab) are expressed by 



61 = 1 



Wh 0) 



o^-co% 



, ^2=- 



^,^3=1- 



(X) (x)^ — 0)t 



(3) 



where the sign of coh is determined as being identical to the sign of the d-c magnetic field. The 
plasma parameters ej, €2, and €3 in the surrounding region are given by (3) only by replacing ojp 

by w;. 

When the z-dependence of the electromagnetic fields is assumed to take the form of e~^^^^^ 
(h is the relative axial propagation constant), the field components are obtained by a combination 
of two scalar functions which satisfy the wave equation [Kales, 1953]: 
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(4) 



where Si and 52 correspond to relative transverse propagation constants (real part is an attenuation 
constant) of a left-handed plane wave and/or a right-handed plane wave (any combination depend- 
ing on the value, A^) for a Cartesian coordinate, and these two values are mutually coupled through 
the value h^ :sl and si are the two solutions of the quadratic equation with respect to s^ given by 



€1/ Ci 



+ f^[{h'-iO'-€l]=0. (5) 



The relative propagation constants 5J, ^2 in the surrounding region also satisfy the equation of the 
same form with (5) in which €1, €2, and 63 are replaced by ej, e!^, and e'j. From this equation we 
obtain the following relations: 



s;^+5^^=A^(i+|)+^-6;-6; 



^1 



(6) 



which relate the transverse propagation constants in the surrounding region to the axial propaga- 
tion constant. 

As is well known, the electromagnetic fields can neither be separated TE and TM modes 
nor £■ — symmetry modes and //—symmetry modes in general. For the present problem we obtain 
the electromagnetic fields in the following forms: 



E\ =j\ {As[(s';' + K'^')e-''^^'^^Bs:,{s[^^ K''')e-'^^''^} 
^!r =j ^ {Ae-''^^'-" -h ^e-^^'^"^) 



^ Zq/I 



(7) 






H\ = -^{As [e-^'^f^^-^ + Bsl.e-'^^^''], 



x^ a 



for the upper half of the surrounding region (signs of s[ and 5^ are changed for the lower half, 
X ^ — a), and 

E^^ = c(si^X'-\-^] cosh sikox-\-D(sl-hjV^^j cosh 52^:0^ 
E^y =—jey(C cosh s\k()X + D cosh s-ikox) 



€1 



E^^ =—j-\{Cs\(s\^ K^) sinh 5iA:oX-h 1)52(5^ +A^) sinh 52/r():^} 
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(3) 



//" =7^-(C cosh Sikox-\-D cosh S2kox) 
H'^=^{C{s^-^A^)coshsikoX-\-D{sl-\-A^) cosh S2kox} 
H^^ =-^(Csi sinh Sikox-\- Ds2 sinh szkox), 



€i 

for the slab region, respectively. Equation (8) corresponds to the case where the transverse 
field components are symmetric with respect to the y — z plane, while the axial field components 
are antisymmetric. This equation gives the modes which reduce to a whistler plane wave prop- 
agating in the z-direction when the thickness of the slab decreases to zero, or ojp becomes equal 
to (Op. Another mode whose transverse field components are antisymmetric with respect to the 
y^-z plane are expressed by the same equation as (8) only by interchanging the hyperbohc sine 
functions and hyperbolic cosine functions. 

By assuming the conventional boundary conditions to be valid at the interfaces, x=±a^ we 
obtain the rigorous dispersion equation expressed as follows: 

2 ^^^^ h'iXt -X2) - {s[ + 4)(5? - 5i)(e,;Z,Z2 + €ss[s:,) - €',(s[s', - A'2){(52 + A2)Zi - {si + A2)Z2} 
€161 

-es{s';' + s!/-hs[s!,^X''){{sl-hX')X2-{sl-^A')Xi}=0, (9) 

where 

Xi,2 = si,2 1 (5i,2A:oa), 

and where Zi, 2 takes a hyperbolic tangent function for the case of symmetric transverse fields 
and a hyperbolic cotangent function for the case of antisymmetric transverse fields. 
For later use the ratio of B to ^ is given below: 

S;(sf-5|) + (52 + A2)J2-(5i+A2)Z,+^(5;2 + A'2)(Z,-J2) 
— ^'^ g(s^-s;)fcoa. (10) 



4(sf - si) + (si + A2)J2 - (4 + A')Xt + ^ (s!/ + A'^)(Xi -X2) 



It can be easily shown from (9) and (10) that, when the thickness of the slab decreases to zero, 
or a)p approaches to co^, s[ becomes zero, and therefore F and the axial field components vanish. 



4. Approximate Solutions for a Thin Slab 

In this section approximate but exphcit solutions are derived for the slab whose thickness 
is small compared to a free space wavelength, and for the hyperbolic tangent mode (Zi,2 in (9) 
takes the tanh function) which reduces to the whistler plane wave when koa decreases to zero. 
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If koa is set equal to zero we have s[ (or 52)^=0 from (9) for the tanh function, and therefore 
from (6) we obtain a whistler plane wave solution: 

^, (11) 

s^^ = ±-f(e;±e.He;,). 

By making use of the above relation we can solve the dispersion equation approximately as a per- 
turbation of a whistler plane wave propagating in the 2-direction. Let koa be sufficiently small 
so that the relation, 5i, 2koa <^ 1 holds however large 5i, 2 is. Neglecting the terms of (koa)^ and the 
higher orders in (9), we obtain the following approximate dispersion equation: 

€3S[s;^(s[-\-S;^) + koafis[, 5;, 51, 52, h) = 0, (12) 

where 

/-e.;(5;5^-A'2) (52+52+ A2) -63 A2(5;2+5^2^5;5: + A'2)-2^^^'/i2. (13) 

^i€, 

Substitution of 5^ = in (13) leads to a rather compact form: 



/0-^^^{€?-(6;±6;=P€2)n, (14) 



after lengthy algebraic calculations. Thus, s[ is obtained from (12) and (14) as 

5, = 7^ Koa 



^'r^-^Wi- {e[±€',^e2)^} koa^ (15) 



61 (6;± €;+€;) 



It should be noted here that (15) is not vahd when |52o| is zero or too small, i.e., €^±€2 + 63 — 
(actually this occurs at co — a)///2). By making use of (15) and the relation that 5 J should satisfy 
(5) for the primed e's, we obtain h^: 



h^ = hl + '-^^f^s[ 



2€. 



and therefore from (6) Sj^ is obtained as 

4' = s;,l + ^{€[^ + e','±e'M + 4)}s['. (17) 

In the preceding equations the upper sign corresponds to the perturbed wave mode of the 
left-handed plane wave and the lower sign corresponds to the one of the right-handed plane wave 
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or the whistler mode when co < cuh- Now we shall restrict ourselves to the whistler type of mode 
and assume that 

o>^,co;^>a^. (18) 

This assumption is considered to be valid in the magnetosphere at most cases. Under the con- 
dition of (18) s[ and s^^ are expressed as follows: 

, {X-X'){X+X'-Y(X -X')} , ^^^ 
s[=- X{Y-2) —koa,Y^2 (19) 

where 

By examining the behavior of the above two constants when the constitutive parameters X, 
X\ and Y are varied, we can know the existence of the trapped or partially trapped waves and 
their properties. Equation (19) shows that s[ is a real number, at least in the first order approxi- 
mation. In order that the waves are trapped along the slab, 5^ must be positive so that the fields 
decrease exponentially in the transverse direction outwards from the slab. From (19) we find 
that there exist three cases for s[ > 0: 

Case (i): X <X' and (j)hi2 < o) < ojh 

Case (ii): X> 3X' and a>///2 < (o < coc (21) 

Case (iii): X' < X < 3X' and ojc < o) < (Oh/2, 



where 



From (20) we find that 



X — X' 

o)c = ajH ^_^^, ' (22) 



,^ > 0, a)Hl2 < C0< 0)H, 
^' I <0, 0< (0<WHI2. ^^^^ 



'2 



It follows that 52 is a real number, i.e., S2=± 1^2 1 for Cases (i) and (ii), and is an imaginary number, 
i.e., 52==±7|52| for Case (iii). It is possible to show from (5) that 52 is a positive real number 
when h^ is given by (16) for both Case (i) and Case (ii). Therefore we find that these two cases 
give the completely trapped surface wave modes whose fields decay exponentially outwards from 
the slab. In Case (i) the plasma density of the slab is less than that of the surrounding region. 
And Case (ii) corresponds to the case where the plasma density of the slab is greater than three 
times of that of the surrounding region. The upper cutoff frequency given by (22) of this mode 
increases with increasing enhancement of the plasma density of the slab. 

In Case (iii) where the electron density of the slab is between one and three times of that of 
the surrounding region, an imaginary value for 53 violates the principle of energy conservation 
since it represents a whistler plane wave outgoing from or incoming to the slab whereas h^ is 
a positive real number representing a nonattenuating wave. In fact a pure imaginary value for 
52 with the other propagation constants retained as real numbers can not satisfy the rigorous 
dispersion equation (9). Thus (9) must be solved more rigorously as a complex dispersion equation 
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in order to be consistent with the energy conservation. It is believed however that the negative 
values of ^2^ represented by (17) or (20) are valid in the first order approximation for the phase 
constant, leading to the leaky surface wave which is partially trapped along the slab and leaks 
part of its energy continuously through propagation. The attenuation constant due to the leakage 
of the power and hence the imaginary part of s'^ is found to be in the order of (/coa)"^ from the 
following discussion. 

Without recourse to solving the exact dispersion equation for complex values of //^, we may 
be able to obtain an approximate attenuation constant of the mode for Case (iii) by power con- 
sideration, that is, by obtaining the power transmitted along the slab and the power leaked from 
the slab. When k^a is sufficiently small compared with unity the transmitted power is substantially 
contributed by the field components associated with s[ in the surrounding region. Thus we 
obtain the approximate transmitted power per unit length in the y-direction {A = 1) as 

P, = 2 Re f " t^. X H* -£, X H*]dx = 2 -^- (24) 

jo 5j/CoZo 

The transmitted power increases in inverse proportion to s[ or a for a constant A. The leaked 

power per unit axial length can be obtained by using the fields at a far distance from the slab 

where the field component associated with s\ diminishes. Thus we obtain the approximate 

leaked power per unit area of the slab: 

P, = 2Ke[£,X//*-E.X//*].^^ ^-2;|^(6i-6: + 6;)|.s^||^|^ (25) 



in which equation the sign of 6^ (•^2'' ^ ^^ '^ determined as ^o = — y'l-s.^l so that only outgoing (energy) 
wave from the slab is permitted from the physical reality: the phase is directed towards the slab 
conversely with the power in this case [Smith, 1960J. Thus, the normalized attenuation constant 
a with respect to Ao is obtained as 

<^ = -fK''s[\r\HY-2)\/Y{Y-2), (26) 

where F and s[ are given by (10) and (19) respectively, and both are proportional to koa (<^ 1), it 
follows therefore that a is proportional to (Aoa)'^ We find also by inspecting (26), (10) and (19) 
that a increases very rapidly when Y decreases towards 2. 

5. Exact Solutions 

The rigorous dispersion equation given by (9) has been numerically solved by using an elec- 
tronic computer for Case (i): o)p < co'p and co///2 < co < co//, where the completely trapped surface 
wave modes have been found to exist in the preceding section. In the same frequency region 
the completely trapped nonattenuating surface waves whose transverse fields are antisymmetric 
with respect to y—z plane (that is, coth function is taken as Zi,2 in (9)) are also found to exist, if 
the plasma frequency of the slab is below a certain critical value o)^ which is less than o)' depending 
on the frequency. This mode shows the dispersion property quite similar to the tanh modes. 
The calculated values of h^ are shown in figure 3 for the tanh function and in figure 4 for the 
coth function in the next section, where the approximate solutions obtained in the preceding 
section are compared with the exact solutions. 

6. Numerical Examples and Discussions 

Figure 2 shows the approximate phase constant for Case (i). The parameter is the thickness 
of the slab in terms of the plasma wavelength of the surrounding region. In this figure the nu- 
merically obtained exact values (only for a pure real h) are plotted by crosses and compared with 
the approximate values. The excellent agreement is found between them. Due to the vanishing 
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Figure 2. Approximate relative phase constant of the 
trapped wave mode (Case (i)), Wh/o)' =0.3, v = XIX' 
= 0.81, Xp is a free space wavelength corresponding 
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of the denominator of the expression for h^ at co = aj///2, the approximation is invahd near 1/7=0.5, 
where the exact h becomes complex. 

In figure 3 the exact solutions obtained numerically for the same mode above (tanh mode) 
are plotted. We find from this figure that the moderate inhomogeneity due to the slab gives rise 
to considerable deviation of the phase constant from that of the plane wave. This is explained 
by the fact that most of the transmitted energy is confined in the slab with a plasma density less 
than that of the surrounding region. 

Figure 4 shows the phase constant for the coth mode obtained by using a computer. The 
dispersion curves are quite similar to those of tanh mode in figure 3, but show less deviation from 
the curve for the plane wave than in the case of tanh mode. This indicates less confinement of 
energy than the tanh mode as seen from the field configurations of these two modes. 

Figure 5 shows the phase constant for Case (iii): the partially trapped (leaky) wave mode. 
The lower cutoff frequency coc given by (22) exists in this mode. Below this frequency s[ is nega- 
tive and the slab does not act as a guiding medium. But this does not necessarily imply that 
the propagation is forbidden: figure 7 shows that the attenuation decreases with decreasing 
frequency. The approximation is again invalid near the frequency, co///2. We find that the 
deviation of the phase constant from the values for a plane wave is rather small compared to that 
in Case (i) (figure 2). This coincides with the fact that the wave is partially or loosely trapped to 
the slab. 

Figure 6 illustrates the transverse field distributions of the partially trapped surface wave. 
The field amplitudes decreases almost exponentially and very gradually near the slab with ripples 
which are caused by the interference between the trapped wave component and the leaked wave 
component. The amphtude ratio of these two wave components is constant as determined by 
(10) to satisfy the boundary conditions. It follows therefore that the trapped wave component 
is continuously transformed into the leaky wave component through the mutual power between 
them, and this gives rise to an attenuation of the transmitted power. This mode is improper 
in the sense that the fields increase to infinity theoretically when x-^±<x) if the surrounding 
plasma is lossless. This is a common nature of leaky waves. From figure 6 it is found that 
the waves are more strongly trapped with increasing frequency and thickness of the slab. 

Finally figure 7 shows the attenuation of the partially trapped wave due to the leakage in 
nepers per one free space wavelength. We find that the attenuation increases with increasing 
frequency, thickness of the slab and enhancement of the plasma density of the slab. The thicker 
slab causes the stronger couping between the forementioned two wave components and results 
in greater attenuation. As an example, the attenuation is about 3.0 dB when /= 5 kc/s, /// = 50 
kc/s,y^ = 500 kc/s, p = X/X'=l,03, the slab of 0.34 km in width, and the propagation path of 10^ km 
are assumed. 
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Figure 3. Exact relative phase constant of the trapped 
wave mode (tanh mode of Case (i)), Wh/co' = 0.3. 
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Figure 5. Approximate relative phase constant of the 
partially trapped wave mode (Case (Hi)), o)hI(o^ 
= 0.33, v=l.23. 
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Figure 4. Exact relative phase constant of the trapped 
wave mode (coth mode of Case (i)), oju/cxj' = 0.3. 



Figure 6. Field distributions of the partially trapped 
wave mode in the transverse direction, (n^joi' =0.33, 
v = 1.23, k' = 27r/k'. "" 



Figure 7. Attenuation of the partially trapped wave 
mode due to the leakage of power in nepers per one free 

space wavelength, cji^loj' = 0.1, ; v = l.03, ; 

v=l.L 
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7. Conclusion 

A full wave theory has been applied to the whistler radio wave propagation along a magneto- 
plasma slab with an enhanced or depressed plasma density which is imbedded in an infinite mag- 
netoplasma. The approximate but explicit solutions for the thin slab have been presented, and 
part of these approximate solutions have been confirmed by the numerically obtained exact solu- 
tions. The existence of the trapped or partially trapped waves and their dispersion properties, 
the fields, the attenuation due to the leakage and others have been discussed. The principal 
results obtained are as follows: 

(i) The propagation of the completely trapped surface wave modes is possible for the slab 
of the depressed ionization in the frequency region of a)///2 < co < coh- 

(ii) The same type of surface waves is propagated also along the slab whose plasma density 
is greater than three times of that of the surrounding plasma, in the frequency region of a)///2 

< 0) < OJc. 

(in) In the frequency region of coc < w < w///2 only a partially trapped (leaky) surface wave 
mode is propagated along the slab of the enhanced ionization. The cutoff frequency ooc depends 
on the degree of the enhancement of ionization. 

(iv) The attenuation due to the leakage of the partially trapped wave increases very rapidly 
with increasing frequency, thickness, and enhancement of ionization of the slab. 

The above results have been derived for the sharply bounded slab of homogeneously en- 
hanced or depressed ionization because of relative simplicity of the analysis. It should be noted 
therefore that the quantitative conclusions in the above can not be strictly applied to an inhomo- 
geneous slab with continuous boundaries, particularly when the thickness of the slab becomes 
appreciable in terms of a free space wavelength, and thus a different approach will be needed 
for such models [Wait, 1962]. 



The author expresses his gratitude to Y. Mushiake of the Tohoku University for his continuing 
interest and advice in the course of this work. 
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